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Gonzalez  RR,  Cheuvront  SN,  Ely  BR,  Moran  DS,  Hadid 
A,  Endrusick  TL,  Sawka  MN.  Sweat  rate  prediction  equations 
for  outdoor  exercise  with  transient  solar  radiation.  J  Appl  Phys¬ 
iol  112:  1300-1310,  2012.  First  published  January  12,  2012; 
doi:  10.1 152/japplphysiol. 01056. 201 1. — We  investigated  the  valid¬ 
ity  of  employing  a  fuzzy  piecewise  prediction  equation  (PW)  [Gon¬ 
zalez  et  al.  J  Appl  Physiol  107:  379-388,  2009]  defined  by  sweat  rate 
(msw,  g-m^-h”1)  =  147  +  1.527-(Ereq)  -  0.87  -(Emax),  which 
integrates  evaporation  required  (Ereq)  and  the  maximum  evaporative 
capacity  of  the  environment  (Emax).  Heat  exchange  and  physiological 
responses  were  determined  throughout  the  trials.  Environmental  con¬ 
ditions  were  ambient  temperature  (Ta)  =  16-26°C,  relative  humidity 
(RH)  =  51-55%,  and  wind  speed  (V)  =  0.5-1. 5  m/s.  Volunteers  wore 
military  fatigues  [clothing  evaporative  potential  (im/clo)  =  0.33]  and 
carried  loads  (15-31  kg)  while  marching  14-37  km  over  variable 
terrains  either  at  night  ( N  =  77,  trials  1-5)  or  night  with  increasing 
daylight  (N  =  33,  trials  6  and  7).  PW  was  modified  (Pw,sol)  for 
transient  solar  radiation  (Rsoi,  W)  determined  from  measured  solar 
loads  and  verified  in  trials  6  and  7.  PW  provided  a  valid  msw 
prediction  during  night  trials  (7-5)  matching  previous  laboratory 
values  and  verified  by  bootstrap  correlation  (rbS  of  0.81,  SE  ±  0.014, 
SEE  =  ±  69.2  g  m~2  h_1).  For  trials  6  and  7,  Ereq  and  Emax 
components  included  Rsoi  applying  a  modified  equation  Aw, sol,  in 
which  msw  =  147  +  1.527 -(Ereq.soi)  —  0.87 -(Emax).  Linear  prediction 
of  msw  =  0.72 -Pw, sol  +  135  (N  =  33)  was  validated  (R2  =  0.92; 
SEE  =  ±33.8  g-m_2-h_1)  with  PW  [3-coefficients  unaltered  during 
field  marches  between  16°C  and  26°C  Ta  for  msw  £  700  g-m_2-h— k 
PW  was  additionally  derived  for  cool  laboratory/night  conditions 
(Ta  <  20°C)  in  which  Ereq  is  low  but  Emax  is  high,  as:  PW,cool 
(gm~2h_1)  =  350  +  1.527-Ereq  —  0.87-Emax.  These  sweat  pre¬ 
diction  equations  allow  valid  tools  for  civilian,  sports,  and  military 
medicine  communities  to  predict  water  needs  during  a  variety  of  heat 
stress/exercise  conditions. 

thermoregulation;  modeling;  load  carriage;  environmental  indexes; 
fluid  replacement 


numerous  public  health  (21,  38),  military  medicine  (6-11,  23, 
28,  29)  and  sports  medicine  (29,  39)  situations  exist  where  it  is 
important  to  estimate  water  needs.  When  performing  physical 
exercise  in  warm-hot  environments,  sweat  loss  is  the  “critical” 
factor  to  calculate  water  requirements  (18,  21,  23,  38).  We 
recently  developed  from  laboratory  experiments  a  fuzzy  piece- 
wise  equation  (PW)  that  predicts  measured  sweat  losses  within 
a  standard  error  estimate  (SEE)  of  ±137  ml/h  (17).  The 
equation  was  also  validated  against  several  laboratory  (indoor) 
studies  with  soldiers  wearing  military  clothing-equipment,  as 


Address  for  reprint  requests  and  other  correspondence:  S.  N.  Cheuvront, 
U.S.  Army  Research  Institute  of  Environmental  Medicine,  Kansas  St.,  Bldg.  42, 
Natick,  MA  01760  (e-mail:  samuel.n.cheuvront@us.army.mil). 


well  as  one  field  study  (outdoors)  in  which  soldiers  wore 
chemical  protective  clothing  while  exposed  to  constant  and  low 
solar  loads  and  low-intensity  and  high-intensity  metabolic  rates 
(M,  W/m2).  Therefore,  both  the  initial  model  development  and 
the  validation  studies  represented  conditions  impacted  by  con¬ 
stant  mild-to-moderate  solar  radiation. 

Sweating  rate  (msw,  g  •  m  2  ■  h  1 )  is  calculated  from  PW  as: 
147  +  1.527-(Ereq)  —  0.87-(Emax)  (17).  Prediction  accuracy  is 
improved  by  50-60%  over  legacy  equations  (41)  and  PW  has 
wide  applicability  to  a  broader  range  of  environmental  temper¬ 
atures,  metabolic  rates,  and  modern  military  clothing  equip¬ 
ment  configurations.  However,  it  is  unknown  if  the  PW  equa¬ 
tion  provides  an  improvement  over  other  conventional  equa¬ 
tions  that  predict  sweat  loss  with  solar  load  (42).  Furthermore, 
PW  has  not  been  validated  over  extensive  military  operational 
missions  (outdoors)  that  encompass  heat  loads  from  transient 
solar  radiation.  Previous  studies  have  shown  that  sweat  rate  is 
altered  by  solar  heat  flux  (14,  30)  and  the  whole  body  sweat 
rate  is  elevated  proportional  to  the  magnitude  of  thermal  and 
exercise  load.  Accurate  measurements  of  required  evaporative 
heat  loss  (Eieq)  can  be  predicted  by  analyzing  key  heat  ex¬ 
change  parameters  of  the  heat  balance  equation  (5,  15,  20,  24). 
The  impact  of  solar  radiation  in  modifying  heat  exchange,  and 
thus  sweat  rate  prediction,  is  easily  determined  empirically  (15, 
16,  24,  25,  30).  The  most  strenuous  military  operations  are 
often  conducted  in  early  morning  hours  that  are  subject  to 
transient  solar  radiation  coupled  with  cooler  environmental 
temperatures.  A  requirement  exists  to  carefully  measure  sweat 
rate  during  outdoor  military  activities  with  and  without  solar 
radiation  to  determine  the  validity  of  the  PW  equation  for  these 
conditions  and  ascertain  if  PW  should  be  modified  from  ap¬ 
propriate  empirical  measurements. 

The  purposes  of  this  study  were  7)  to  validate  the  original 
PW  sweat  prediction  equation,  derived  from  laboratory  studies, 
for  implementation  outdoors  at  night  and  with  transient  solar 
circumstances  in  soldiers  performing  military  activities  and 
wearing  military  clothing-equipment;  and  2)  to  decide  whether 
a  modified  PW  equation  should  be  developed  de  novo,  and  if 
so,  then  determine  if  it  provides  a  reliable  prediction  of  sweat 
rate  during  prolonged  work  with  transient  solar  radiation  con¬ 
ditions.  Our  hypothesis  was  that  regression  coefficients  in  the 
fuzzy  piecewise  algorithm  deriving  the  PW  (15,  44)  would  be 
robust  enough  to  permit  valid  msw  predictions  during  outdoor 
night  marches,  but  would  require  adjustment  (Pw,sol)  for 
transient  solar  radiation  exposures  and  possibly  cooler  condi¬ 
tions.  Consequently,  an  additional  hypothesis  was  that  the  PW 
equation,  adapted  with  required  solar  load  factors,  will  provide 
a  reliable  sweat  loss  (water  needs)  prediction  equation  that  is 
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Table  1.  Description  of  conditions  for  the  night  marches  ( trials  1-5) 

1301 

Trial  N 

Section 

Time, 

min 

Distance,  km 

Terrain  Type,  %  of  Section 

Terrain 
Factor  (t)) 

1  and  2  25 

i 

60 

6 

Base  course,  40%;  sand,  60% 

1.4 

2 

50 

5 

Hard  road,  40%;  plowed  field,  10%;  hard  road,  50% 

1.4 

3 

50 

5 

Hard  road,  40%;  plowed  field,  10%;  hard  road,  50% 

1.4 

4 

68 

5 

Base  course,  40%;  sand,  60% 

1.4 

5 

30 

3 

Base  course,  40%;  sand,  60% 

1.4 

3  12 

1 

60 

6 

Hard  road,  40%;  wet  hard  sand,  60% 

1.5 

2 

50 

5 

Hard  road,  40%;  plowed  field,  10%;  hard  road,  50% 

1.5 

3 

65 

6.5 

Wet  hard  sand,  80%;  dry  sand,  20% 

1.5 

4 

65 

6.5 

As  above 

1.5 

5 

55 

5.5 

Same  as  section  2 

1.5 

6 

75 

7.5 

Dry  sand,  80%;  1  km  on  hard  road,  20% 

1.5 

4  and  5  40 

1 

55 

6 

Trial  4:  —7.3  km/h  (%gradient  —0.7)  descent  on  gravel,  15  min  rest;  trial  5:  —9.6  km/h  (%gradient 
—0.8),  15  rest 

1.3 

2 

55 

6 

Trial  4:  10  km/h  ascent  on  gravel;  17  min  rest;  trial  5:  7.8  km/h  (%  gradient  0.7)  ascent;  17  min  rest. 

1.3 

3 

28 

2 

Trial  4:  —16.1  km/h  (%gradient  —0.8)  descent  on  gravel,  completion;  trial  5:  1.7  km/h  ascent 
(%gradient  0.2),  march  completion. 

1.3 

N  is  number  of  volunteers  completing  each  trial.  Section  is  the  type  of  road  during  each  march,  %  time,  and  various  ascent/descent  (km/h;  %gradient)  marches 
of  trials  4  and  5.  Terrain  factor  (iq  )  is  unit  coefficient  of  Pandolf  equation  prediction  (31,  32),  as  modified  by  Santee  et  al.  (35)  for  specific  descent  marches. 


functional  for  a  wide  range  of  indoor  and  outdoor  operational 
scenarios  with  improved  prediction  over  legacy  equations  (41, 
42)  and  other  modeling  approaches  (3,  24,  46)  that  account  for 
solar  radiation.  As  in  the  previous  study  (17),  we  applied 
Tseng’s  and  coworkers  (44)  fuzzy  piecewise  analysis  to  predict 
sweat  losses  measured  from  field  data  assuming  that  msw 
(lumped  efferent  eccrine  output)  has  both  linear  and  nonlinear 
characteristics  as  a  function  of  core  temperature  (Tcore)  and 
skin  temperature  (Tsk)  inputs  (15,  24,  37).  Nonthermal  factors 
(43)  that  may  affect  the  general  gain  of  the  efferent  thermal 
drive  and  sweating  threshold  were  not  accounted  for  in  the 
analysis.  The  sweat  loss  equations  (indoor  and  transient  solar 
load)  in  the  present  study  should  be  considered  as  lumped 
derivatives  and  include  central  nervous  system  sudomotor 
drive  principally  from  thermal  factors  with  nonthermal  factors 
having  minimal  effect  (37). 

METHODS 

One-hundred  ten  male  soldier  volunteers  participated  in  numerous 
outdoor  field  trials.  The  appropriate  institutional  human  use  review 
boards  approved  the  protocol,  and  all  volunteers  provided  written 
informed  consent.  Investigators  adhered  to  policies  for  protection  of 
human  subjects  as  prescribed  in  US  Army  Regulations  70-25  and  US 
Army  Medical  Research  and  Materiel  Command  Regulation  70-25. 
The  research  was  conducted  in  adherence  with  the  provisions  of  Code 
45  of  Federal  Regulation  part  46. 

Experimental  studies.  Experimental  data  were  collected  by  staff  of 
the  Heller  Institute  of  Medical  Research.  Seven  outdoor  marches  were 
performed  by  active  military  units  and  incorporated  into  their  training, 


with  accommodations  made  for  careful  data  collection.  The  military 
units  maintained  normal  leadership,  training  hours,  clothing-equip¬ 
ment,  and  work  rates  (load  and  speed)  applicable  to  actual  missions. 
Exercise  level  per  se,  state  of  hydration,  and  sleep  state  in  the  subjects 
all  were  controlled.  Subjects  were  well-hydrated  at  the  beginning  and 
during  rest,  and  none  of  the  subjects  was  sleep  deprived. 

The  outdoor  marches  were  conducted  near  field  stations  exhibiting 
diverse  terrains  (desert,  Mediterranean  coastal,  and  various  foothill 
areas).  Tables  1  and  2  describe  the  number  of  subjects  and  specific 
conditions  for  the  night  ( trials  1-5 )  and  transient  solar  radiation  ( trials 
6  and  7)  outdoor  marches,  respectively.  Data  were  collected  on  77 
volunteers  who  marched  at  night  ( trials  1-5),  and  some  of  these 
individuals  participated  in  multiple  trials.  Thirty  three  supplementary 
volunteers  marched  during  the  transient  solar  radiation  conditions 
(trials  6  and  7).  Table  3  shows  many  of  the  variables  and  equations 
with  essential  modifications  used  in  the  study. 

Physiological  measures.  Oxygen  uptake  was  measured  from  a 
cohort  of  the  volunteers.  These  data  were  then  used  to  estimate 
metabolic  rate  of  other  individuals  in  each  group.  Metabolic  rate  was 
calculated  from  a  90-s  sample  of  expired  air  collected  using  a  K4B 
(COSMED,  Italy)  metabolic  measurement  device.  Marching  pace  and 
grade  were  estimated  using  a  Garmin  GPS  navigator.  Oxygen  uptake 
calculations  were  estimated  based  on  weight,  load,  walking  speed  + 
grade,  and  terrain  factor  using  Pandolf  et  al.  (32)  equations  and  Santee 
et  al.  (35)  modifications  for  some  downhill  trials  to  estimate  net  heat 
production.  In  some  cases,  a  median  terrain  factor  coefficient  was 
chosen  since  portions  of  some  marches  occurred  on  gravel,  sand,  and 
paved  road  (Tables  1  and  2).  The  oxygen  uptake  calculation  was  for 
the  exercise  portion  only,  as  rest  periods  were  estimated  at  a  metabolic 
rate  of  120  W  (average  of  63  ±  3  W/m2)(15,  24).  A  standard 
respiratory  exchange  ratio  for  a  mixed  diet  was  estimated  for  consis- 


Table  2.  Description  of  conditions  for  the  transient  solar  load  marches  (trials  6  and  7) 


Trial 

N 

Section 

Time,  min 

Distance,  km 

Terrain  Type/%  of  Section 

Terrain  Factor  (t)) 

6 

20 

i 

65 

6.0 

Gravel,  100%;  18-min  rest  period 

1.30 

2 

48 

4.4 

Gravel,  100%;  17-min  rest  period 

1.30 

3 

62 

6.0 

Gravel,  100% 

1.30 

7 

13 

1 

58 

5.8 

Paved  road,  43%;  gravel,  57% 

1.25 

2 

71 

6.7 

Gravel,  100% 

1.25 

3 

62 

5.8 

Paved  road,  43%;  gravel,  57% 

1.25 

Abbreviations  are  same  as  Table  1.  Trial  6  began  at  0330  and  ended  at  0700;  average  global  solar  load  (gSL)  monitored  from  2  stations  =  63  W/m2.  Trial 
7  began  at  1420  and  ended  at  1800;  average  gSL  =  265  W/m2. 
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Nomenclature  and  description  of  key  variables  and  equations  implemented 

Variable 

Units 

Description 

Formula 

Refs. 

Modification  in  Present  Study 

PW 

Sweat  rate  prediction  by  fuzzy 
piecewise  (PW)  change  point 
regression 

147  +  (1.527 -Ereq)  - 
(0.87 -Emax) 

17,  44 

Calculation  by  bootstrap  correlation  of 
original  and  night  field  data: 
equation  remains  valid  N  >  500 
cases;  outliers  do  not  statistically 
change  comparison;  SEE  within 
±69  g*m_2*h_1 

M 

W/m2 

Metabolic  heat  production 

Calculated  from  measured 
oxygen  uptake  (V02) 
and  Pandolf  equation 

15,  32,  35 

Pandolf  equation  adjusted  by  Santee 
coefficients 

P  w,sol 

g*m_2*h_1 

Sweat  rate  prediction  for  transient 
solar  load 

147  +  (1 .527  *  Ereq.sol)  “ 

(0.87 -Emax) 

This  study, 
Fig.  3 

Solved  for  Rsoi  per  march  time  using 
transient  solar  field  data;  sweat  rate 
becomes  predictable  as  OLS 
equation:  0.72*Pw,sol  +  135, 

R2  =  0.92 

Emax 

g-m^-h-1 

Maximum  evaporation  capacity  of 
a  given  environment 

Emax  =  U)he(Ps,sk  —  Pa) 

15 

Adjusted  by  skin  wettedness  (co).The 
fraction  of  the  total  body  area  (Ad) 
covered  by  sweat  [wetted  area 
(Aw)],  i.e.,  Aw/Ad.  and  evaporative 
heat  transfer  coefficient,  he 

Ereq 

g*m~2*h_1 

Required  evaporation  from  heat 
balance  includes  total  skin 
evaporation,  all  respired,  non 
eccrine,  and  metabolic  water 
losses 

Ereq  =  M  -  (WeXt)  - 

DRY  -  S 

15 

Valid  for  indoor  lab  and  night  field 
data;  for  transient  solar  field  study, 
Ereq.sol  first  solved  for  Rsoi  as  input 
to  heat  balance  equation  to 
determine  Pw,sol  . 

R  +  C 

W/m2 

Dry  or  sensible  heat  loss 

[6.45  -(Ta  -  TSk)]/IT 

4.  15 

Eres 

W/m2 

Respiratory  heat  loss 

(Eres  "t” Cres)  =  Ad'M[0.0014 

15 

Modified  by  ambient  temperature  and 

Cres 

W/m2 

Convective  heat  loss 

(34  -  Ta)  +  0.0023- 
(44  -  Pa)l 

water  vapor  pressure,  maximum  Ta 
and  Pa,  evaporative  potential  (Emax), 
respiratory  exchange 
ratio  (R) 

Edif 

W/m2 

Skin  diffusion 

Edif  =  0.05 -Emax 

15 

Hires 

g/min 

Metabolic  heat  loss 

Hires  =  Vo2-(R-0.53) 

17 

Pa 

kPa  or  Torr 

Ambient  water  vapor  pressure 

(RH/100)-exp[  18.6686  - 
(4030.183/Ta  + 

235)];Antoine  equation 

15 

hc 

Convective  heat  transfer 
coefficient 

he  =1.2-[(M  -  50)- 
(Pb/760)]039 

15 

Modified  by  walking,  free  and  forced 
convection,  barometric  pressure 

he 

W-m~2-°C_1 

Evaporative  heat  transfer 
coefficient 

2.2*  hc  or  alternatively,  with 
clothing 

(2.2  *im)/(IT*  0.155) 

(Pb),  Lewis  relation  at  sea  level 
(2.2°C/Torr) 

hr 

W-m_2"°C-1 

Radiative  heat  transfer  equation 

hr  =  5.67- 10~8-0.97-0.77- 
[(Tsk  +273)4  -  (MRT  + 
273)4]/(Tsk  -  MRT) 

15 

Unmodified 

It 

Fcle 

clo 

ND 

Total  clothing  insulation,  sum  of 
air  (Ia)  and  fabric  (Ici)  layer 
insulation 

Effective  clothing  factor  for  heat 
exchange  from  articulated 
manikin 

(la  +  Icl) 

4.  15 

clo  =  0.155°C*m2/W  or  thermal 
conductance  of  6.45  W*m2/K 

Ps,sk 

kPa  or  Torr 

Skin  saturation  vapor  pressure 

Antoine  equation  solved  for  TSk 

15 

Rci 

m2-  K/W 

Intrinsic  (fabric)  thermal 
resistance 

Rci  =  Rt  -  l/[(hc+hr)]-fcl 

1.  15 

Reciprocal  of  thermal  conductance 

Rt 

m2  ■  K/W 

Total  clothing  resistance 

0.1 55 -It 

1.  15 

Modified  by  wind  and  body  motion 

VairOr  V 

m/s 

Ambient  air  movement  (velocity 
on  person) 

Input  variable 

7,  15,  17 

Direct  observations 

Veff 

m/s 

Effective  air  motion  produced  by 
wind  speed  and  M 

Vail  +  (0.004)  (  M-Ad  -  105) 

24,  32 

Unmodified 

CT 

Stefan-Boltzman  constant 

5.67  X  10~8W-m~2-K~4) 

1,  15 

Unmodified 

Wext 

W/m2 

External  work  based  on  activity 

Pandolf  equations 

32 

Santee  (35)  modifications 

Ap 

m2 

Projected  surface  area  exposed  to 
direct  beam  solar  flux 

Fanger,  ASHRAE  equations 

1.  13 

Upright  walking 

feff 

% 

Fraction  of  body  surface  exposed 
to  solar  load  at  a  given  time 
(Ar/AD) 

Fanger,  ASHRAE  equations 

1,  13 

Upright  walking 

ciLW 

W 

Long  wave  absorptivity 

Breckenridge  and  Goldman 
equations 

4,  16 

Upright  walking 

aSW 

W 

Short  wave  absorptivity 

Breckenridge  and  Goldman 
equations 

4,  16 

Upright  walking 

Continued 
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Variable 

Units 

Description 

Formula 

Refs. 

Modification  in  Present  Study 

Rsol 

w 

Matthew  solar  radiant  heat  flow 

-0.0003-  (gSL)2  + 

25 

Adjusted  by  Ad  and  Fci,e  for  inclusion 

0.681 -(gSL)  +  3.136 

into  heat  balance  to  determine 

Ereq.sol  (W/m2) 

ERF 

W/m2 

Effective  radiant  field 

hr  (MRT  -  Ta);@ 

15 

Interchangeable  variables 

[(0.835 -RsoO/Ad] 

gSL 

W/m2 

Global  solar  load 

Direct  weather  station  data; 

25 

Direct  measurement  from  weather 

pyranometer  values 

station  data 

Ith 

Total  solar  irradiance  of  a 

Ro-sinp  +  Rdif 

2,  13 

Meteorological  variables 

horizontal  object 

Me.th 

Radiant  flux  emitted  as  thermal 

Shapiro  et  al.,  outdoor 

42 

Unmodified 

radiation  by  a  surface. 

equations 

MRT 

°c 

Mean  radiant  temperature 

Ta  +  ERF/hr 

14,  15 

Unmodified 

Tsk 

°c 

Mean  weighted  skin  temperature 

Mitchell  equation 

27 

3-Site  direct  measurement 

Vair 

Direct  measurement 

Rd 

w 

Direct  beam  solar  load 

Arens  et  al.,  Breckenridge,  and 

2,  4,  13, 

Direct  or  weather  station  values 

Rdif 

w 

Diffuse  irradiance 

Fanger  equations 

Rref 

w 

Reflected  irradiance 

Ad 

m2 

Body  surface  area  (BSA)  from 

Ad  =  0.202  ■  W0-425  •  H0-725 

1,  15 

Observed  values 

DuBois  equation 

tency  based  on  a  respiratory  exchange  ratio  during  submaximal 
intensity  exercise  (50%  Vo2peak)  ranging  from  0.8  to  0.98  (7,  13). 

Although  not  reported,  heart  rate  and  core  temperature  were  mea¬ 
sured  for  medical  monitoring  purposes.  Heart  rate  was  recorded  using 
a  Polar  heart  rate  monitor.  Core  temperature  data  were  recorded  using 
the  telemetric  pill  technique  (HQ  and  Mini  Mitter).  During  the  rest 
periods,  three-site  skin  temperatures  were  measured  by  an  YSI-409 
(YSI,  Yellow  Springs,  OH)  surface  thermal  probe.  Mean  weighted 
skin  temperature  (TSk)  was  calculated  using  a  three-site  (chest,  thigh, 
and  arm)  surface  area  weighting  coefficient  formula  (27).  In  the 
present  data  set,  clothed  Tsk  was  linearly  correlated  with  Ta  by  the 
equation:  Tsk  =  0.53  Ta  +  20.28  ( R 2  =  0.83;  SEE  =  ±0.54)  over  the 
range  of  Ta  from  16  to  26°C.  Skin  saturation  vapor  pressure  (kPa)  was 
determined  for  each  individual  skin  temperature  value  using  An¬ 
toine’s  equation  (15,  24)  and  the  necessary  evaporative  heat  transfer 
coefficients  (/ie)  utilized  to  determine  Emax  for  each  time  interval  in 
the  respective  trials. 

Total  sweat  losses  were  determined  from  changes  in  body  mass, 
corrected  for  nonsweat  losses,  and  assume  that  sweat  volume  and 
mass  are  equivalent  (1  ml  =  1  g).  Sweating  rate  expressed  as  liters/h 
or  g-m_2-h_1  was  determined  by  time  and  surface  area  weightings, 
as  appropriate.  Detailed  calculations  included  (7): 

msw  =  [Anude  body  mass  ^ 

-  (UV  +  EW  +  NEFL)  +  DV]/time 

where  Anude  mass  is  the  difference  in  nude  body  mass  (g)  pre-  to 
postexercise;  UV  is  urine  volume;  EW  is  excrement  weight,  if  any; 
NEFL  is  noneccrine  fluid  losses,  which  include  respiratory  water 
losses  and  CO2-O2  exchange  (7,  36,  37);  and  DV  is  consumed  drink 
volume.  Body  mass  was  measured  while  the  subject  was  minimally 
clothed  and  fully  clothed  before  and  after  each  march  for  estimation 
of  trapped  sweat. 

Environmental  parameters.  The  environmental  parameters  ambient 
air  temperature  (Ta,  °C),  relative  humidity  (RH,  %)  (transformed  to 
ambient  water  vapor  pressure,  kPa),  and  wind  velocity  (V,  m/s)  were 
collected  every  30  min  using  a  Kestrel  3000  (Caliber  Sales  Engineer¬ 
ing)  and  recorded  manually  for  later  spreadsheet  entry.  Composite 
global  solar  radiation  (gSL,  W/m2)  was  recorded  at  two  field  stations 
in  close  proximity  to  the  marches.  The  effective  solar  load  was 
determined  by  calculation  of  Rsol  and  effective  radiant  field  (ERF)  and 
mean  radiant  temperature  (MRT)  (15,  16,  24,  25).  The  procedure 
employed  to  determine  these  coefficients  for  transient  solar  load  will 
be  described  later. 


Shapiro  et  al.  (42)  originally  estimated  solar  radiation  effects  via 
several  correction  factors  annexed  to  their  sweat  loss  prediction 
equation  from  data  collected  during  steady-state  exposure  to  a  variety 
of  solar  load  intensities.  These  factors  were  used  to  estimate  required 
evaporative  heat  loss  (Ereq.sL)  by  assessing  separate  thermal  radiative 
flux  factors  analyzed  from  the  heat  balance  equation  as: 

Ereq,SL( watts)  =  M  -  Wext  +  [Hc  +  Hr  +  HL] .  (2) 

In  the  above  equation  (42),  M  is  metabolic  heat  production  (15,  32, 
35);  Wext  is  external  work  rate  (32);  Hc  is  convective  heat  transfer, 
calculated  as  6.45-AD(Ta  —  Tsk)/IT;  AD  is  the  Dubois  body  surface 
area  (m2)  (15);  and  It  (4)  is  the  fixed  value  of  total  clothing  insulation 
[originally  evaluated  by  a  static  manikin  that  include  the  air  (Ia)  and 
intrinsic  clothing  layers  (Ici)] .  The  radiative  heat  flux  factors  comprise 
Hr  (W),  which  was  estimated  from  radiative  heat  transfer  = 
[1.5  ■  Ad  -  (SLg°-6)/IT],  where  SLg  incorporates  the  shortwave  radiative 
load;  and  long-wave  emission  (HL,  W)  from  the  total  body  to  the 
environment,  which  was  estimated  as: 

Hl  =  (0.047  -Ad-M^/It  (3) 

where  the  factor  (Me,th)  is  the  radiant  heat  flux  that  includes  the 
combined  theoretical  perfect  black-body  long  wave  emission  from  all 
body  surfaces.  Me,th  is  equal  to  ct,  the  Stefan-Boltzman  radiation 
constant  (5.67  X  10~8  W ■  m~2 ■  °K~4)_multiplied  by  T4,  calculated  as 
the  outer  surface  layer  temperature  (Tci)4  and  transformed  into  °K 
(273.15  +  °C)  (4,  15,  25).  Effect  of  heat  storage  (S)  was  not 
accounted  for  in  the  above  analysis,  and  steady-state  was  assumed  for 
all  the  heat  exchange  variables. 

Several  thermal  radiation  algorithms  were  derived  in  the  present 
study  quantifying  a  more  inclusive  radiative  heat  transfer  analysis 
possible  for  the  transient  state  (1,  2,  4,  16,  25).  These  analyses  were 
applied  for  evaluating  transient  solar  load  on  the  individual  [i.e.,  when 
the  global  solar  radiation  [gSL  =  (Hr  +  HL)]  is  not  constant).  Solar 
radiation  is  generally  measured  as  the  sum  of  the  total  (direct  plus 
diffuse)  radiation  falling  on  a  horizontal  surface  (gSL).  In  meteoro- 
logic  terms  this  is  expressed  as  Ith.  which  varies  with  solar  elevation 
(1,  2,  16).  In  addition,  the  area  of  the  human  body  exposed  to  solar 
radiation,  and  therefore  the  solar  energy  received  by  the  body,  also 
varies  with  solar  elevation.  Generally,  gSL  [or  its  sequel,  ITh,  (2,  4)] 
can  be  plotted  for  a  specific  solar  elevation.  A  key  property  is  the 
projected  area  (Ap)  solved  as  a  function  of  the  solar  altitude  and 
normally  is  about  20-25°  at  the  highest  solar  altitude  and  drops  to 
around  80-90°  (1,  2,  13). 
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The  equivalence  between  Ith,  ERF,  and  MRT  (all  inclusive  radi¬ 
ative  heat  transfer  algorithms)  involves  the  following  assumptions. 
Radiation  (R)  and  convection  (C)  in  the  heat  exchange  are  assumed  to 
be  equal  in  plotting  MRT.  This  heat  exchange  proportionality  occurs 
because  the  (MRT  —  Ta)  is  directly  proportional  to  ERF  (14).  In  the 
present  study,  long-wave  radiation  characterized  by  MRT  and  ERF  is 
treated  as  uniformly  distributed  (4tt)  on  the  individual  during  each 
walking  phase  (15,  16).  There  exist  no  provisions  for  the  directional 
characteristics  of  the  body’s  exposed  surfaces  interacting  with  non¬ 
uniformities  in  the  long-wave  radiant  field  (4).  The  total  body  surface 
area  of  the  assumed  person  impacted  by  the  solar  load  is  initially 
estimated  by  Ad  (15).  The  projected  area  (Ap)  of  an  individual 
exposed  to  direct  beam  solar  load  while  walking  is  calculated  for  our 
purposes  as  AP/AD  =  0.23  (1,  2,  4,  13-16). 

The  effective  radiant  field  (ERF)  (1,  14,  15)  is  one  convenient  and 
precise  measurement  of  the  net  radiant  heat  flux  to  or  from  the  human 
body.  ERF  incorporates  the  summed  long-wave  radiation  energy 
received  by  the  body  in  a  clothed  state  when  surrounding  surface 
temperatures  are  different  from  the  air  temperature.  Additionally,  the 
surrounding  surface  temperature  may  be  expressed  by  MRT  over  a  4tt 
radiating  area  of  the  body  surface  (14,  15)  particularly  with  sustained 
marching  or  when  surface  temperatures  are  higher  than  extremities  or 
other  microclimate  spaces  surrounding  the  individual.  The  ERF  on  the 
human  body  is  related  to  MRT  by: 

ERF  =  feff  •  hT  (MRT  -  Ta)  (4a) 

in  which  MRT  is  easily  transformed  by 

MRT  =  Ta  +  ERF/(feff  •  hr)  (4b) 

where  feff,  calculated  as  (Ar/Ao),  is  the  fraction  of  the  body  surface 
exposed  to  radiation  from  the  environment;  hT  includes  the  radiative 
heat  transfer  coefficient;  and  Ta  is  the  air  temperature.  Typically,  any 
transient  solar  load  by  ERF  regardless  whether  it  is  direct,  diffuse,  or 
from  the  ground  up  to  the  individual  (W/m2)  divided  by  the  radiative 
coefficients  (m2K/W)  increases  the  effective  temperature  (ET*15) 
above  the  prevailing  ambient  temperature  (°C)  and  intensifies  the  heat 
stress  as  MRT  rises  (15,  16,  25).  Gagge  and  Hardy  (14)  demonstrated 
that  every  interval  of  this  radiative  heat  flux  influences  ERF  and  is 
linearly  correlated  with  thermoregulatory  sweating  rate  assessed  by 
evaporative  heat  loss  from  body  weight  loss  changes. 

For  both  trials  6  and  7,  the  body  surface  areas  of  the  33  walking 
volunteers,  as  sunlight  first  appeared,  were  generally  considered  as 
being  exposed  to  mixed  diffuse,  reflected,  and  long-wave  radiation  in 
which  the  feff  =  0.72  determined  using  Fanger’s  projection  factors 
(13).  The  area  exposed  was  estimated  as  0.72  - AD  =  1.37  m2  over  a 
body  surface  of  1.9  m2  for  the  subjects’  averages.  A  solar  elevation 
(P)  of  45°  was  considered  for  the  walks  in  this  study  (1,  2,  13,  16). 
Diffuse-sky  and  ground-reflected  solar  radiation  (25)  were  assumed  to 
be  uniformly  spread  on  one-half  the  exposed  portion  of  each  of  the 
volunteer’s  BSA  when  clothed  as  feff- 0.95  =  0.684  m2  over  half  of 
the  body  surface. 

Because  ERF  incorporates  the  feff  of  the  receiving  individual,  for 
this  study  the  ERF  integrates  absorbed  energy  from  long-wave  (Hl) 
and  shortwave  (Hs)  sources  originally  found  in  the  Shapiro  (42) 
formulation  estimated  as: 

ERF  ■  aLW  =  [0.72  ■  (0.95/1.9)  ■  (Rdif  +  Rref)  +  (Jfl) 

(0.41/1.9)  RD]  asw 

where,  for  the  walking  individual,  ctlw  is  long-wave  absorptivity,  = 
0.95;  asw  is  short-wave  absorptivity,  =0.67  for  (white)  skin  and  military 
clothing  (for  this  study  IT  =  1.26  clo);  RD  is  direct  beam  solar  radiation 
measured  perpendicular  to  the  beam  (W/m2);  Rdif  is  diffuse  irradiance 
(W/m2)  of  an  upward-facing  horizontal  surface;  and  Rref  is  reflected 
irradiance  (W/m2)  of  a  downward-facing  horizontal  surface. 

Since  RD  =gSL  OR  (Ith  —  Rdif)/sin  p,  assuming  Ith  is  total  solar 
irradiance  of  a  horizontal  surface,  Eq.  5a  can  be  rewritten: 


ERF-0.95  =  [(0.72/2)-(RD  +  Rref)  +  (5fc) 

(0.41/1. 9- sinp)  (lTH  -  Rdif)l  ■  0.67] 

For  the  present  study,  from  Fanger’s  estimates  (13),  Rref  was 
assumed  to  be  0.23  -  Ith  ,  which  is  the  mean  value  of  the  terrain  areas 
during  the  walks  as  solar  gain  (from  a  specific  station  gSF  value) 
accumulates  on  the  individual  during  the  time  periods. 

In  general,  the  average  solar  radiation  load  (as  gSF  or  ITh)  from  the 
various  stations  is  not  always  the  best  approximation  to  quantify 
Ereq.sL  because  the  lumped  quantity  assumes  steady-state  direct  radi¬ 
ant  load  and  each  subject  is  generally  gaining  a  different  net  gSF  per 
time  during  the  walk.  In  this  study,  we  are  interested  in  whole  body 
sweat  response  affected  by  the  gSF  per  each  work/rest  period. 

Solar  parameter  estimation  for  transient  exposures.  Matthew  et  al. 
(25)  extended  the  above  solar  radiation  calculation  approach  by 
developing  an  equation  that  is  applicable  to  the  transient  solar  data  of 
the  present  study  data  analysis.  The  equation  can  be  straightforwardly 
migrated  into  various  heat  strain  models  (24,  31,  46)  that  input 
variables  of  MRT.  ERF,  and  Ta  or  a  combination  of  these  (15,  16,  24). 
Extending  Matthew’s  calculations,  we  developed  a  continuous  func¬ 
tion  that  integrates  both  Hr  and  HL  [from  the  Shapiro  (42)  outdoor 
study]  and  the  Ith  (—  gSF)  concept  discussed  above  to  evaluate 
transient  solar  heat  load  that  each  subject  is  gaining  per  time  during  a 
march.  Our  improved  solution  allowed  us  to  resolve  the  average 
“constant”  solar  radiation  variations.  Rsoi  (W  per  person  in  this  case) 
theoretically  represents  the  best-case  alternative  correcting  variable 
clothing  factors  and  transient  radiant  loads  by  MRT  and  ERF  (4,  15, 
16).  Additionally,  Matthew  and  colleagues  (25)  reported  that  Rsoi  and 
ERF  solutions  are  tightly  correlated  (R2  =  0.93)  over  wide  solar 
domains.  Their  separate  influences  can  therefore  be  interchanged.  As 
such,  MRT  that  is  proportional  to  ERF  becomes  a  convenient  index  as 
an  input  into  various  thermal  models  can  be  easily  solved  from  both 
Rsoi,  ERF,  and  rationally  estimated  from  mean  gSF,  assuming  the 
given  weather  station  uses  standard  meteorological  measuring  instru¬ 
ments. 

In  the  heat  flow  form  expressed  by  Matthew  et  al.  (25), 

Rsol  (watts)  =  -0.0003  •  (gSF)2  +  0.681  •  (gSL)  +  3.136  (6) 

In  the  polynomial  fit  (R2  =  0.93)  of  Eq.  6,  MRT  (°C)  was  derived 
from  measurements  using  a  pyranometer  (to  measure  short-wave 
radiation  fluxes)  and  pyrgeometer  (to  measure  long-wave  radiation 
fluxes).  Normally,  the  value  of  the  solar  load  differences  (MRT  —  Ta) 
is  the  effective  heat  stress  input  parameter  employed  in  various 
thermal  models  (1,  3,  24).  This  parameter  may  be  solved  in  the 
equation  to  estimate  solar  load  on  the  skin  surface  (Hs>sk)  that  is 
currently  implemented  in  the  USARIEM  Heat  Strain  Decision  Aid 
(HSDA)  in  which  Hs,sk  is  computed  as  a  function  of  various  clothing 
factors  (y-clo  and  yi  =  im/clo)  and  Veff  corrected  for  walking 
[“pumping”  factors  (4)  and  load  carriage  variables  (3,  17,  24,  31,  32, 
46)].  In  trials  6  and  7,  MRT  was  calculated  by  solving  for  the  average 
Rsoi  and  its  sequel  (ERF)  at  each  marching  period  and  time  of  day 
from  the  gSF  collected  at  the  various  field  stations.  The  variable  Ith 
will  not  be  considered  further. 

Either  Rsoi  or  ERF  can  be  utilized  in  the  heat  balance  equation  (as 
a  rational  variable  in  the  heat  balance  equation)  by  combining  the 
parameters  Shapiro  (42)  employed  to  evaluate  radiative  short-wave 
(gSF,  W)  and  long-wave  (HLg,  W)  heat  stress  less  radiative  long-wave 
heat  loss  (Hl,  W).  Applying  parti tional  calorimetric  analysis  (15,  16, 
24),  the  heat  balance  equation  (W/m2)  becomes: 

±S  =  (M  -  Wext)  -  Esk  -  (R  +  C)  +  Rsol  (7) 

where  S  is  the  rate  of  heat  storage  less  metabolic  heat  production  (M) 
minus  external  work  rate  (Wext),  minus  radiative  (R)  and  convective 
heat  loss  (C),  and  ESk  =  (E  —  Eres  —  Cres  —  Edif  —  mres).  The  variable 
(E)  incorporates  total  sweat  loss  minus  respiratory  (Eres)  and  convec¬ 
tive  heat  loss  (Cres)  minus  skin  diffusion  (Edif)  less  any  metabolic  heat 
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loss  (mres).  The  Rsoi  in  Eq.  7  is  calculated  as  a  radiant  flux  (per  m2). 
All  avenues  of  heat  exchange  and  respective  heat  transfer  coefficients 
can  be  calculated  as  in  previous  studies  (15,  16,  24). 

The  transient  equation  describing  heat  exchange  between  the  skin 
surface  and  the  radiant  environment  varying  with  time  (Ereq,soi)  (Table 
3)  is  therefore  estimated  by: 

Ereq.soi  (W/m 2)  =  net  heat  flux  -  DRY  + 

effective  solar  gain  ±  heat  storage 

Ereq.soi  (W/m2)  =  Msk  -  hr+c  •  Fcke(Tsk  -  Ta)  + 

[Rsol  ■  Fcl,e]  ±  s 

where  Msk  in  Eq.  8b  now  includes  the  net  heat  flux  from  the  interior 
body  to  the  skin  (M  —  Wext  —  Eres  —  Cres  —  mres).  The  variable  [Fci,e] 
is  the  effective  clothing  factor  determined  from  an  articulated,  mov¬ 
able  manikin  in  separate  analysis  (4,  13,  16,  24). 

In  Eq.  8b,  equivalent  mathematical  substitutions  can  be  used 
assuming  Rsoi  =  ERF  =  hr  •  (Ar/AD)  •  (MRT  —  Ta).  The  variable, 
Ereq.soi,  therefore,  represents  the  solution  of  all  avenues  of  the  heat 
balance  during  transient  solar  loads  over  a  given  time  interval.  The 
value  is  substituted  into  the  PW  equation  and  iterated  into  the  fuzzy 
piecewise  analysis  (44).  An  independent  computer-generated  sensi¬ 
tivity  analysis  showed  that  ERF  (W/m2)  is  thermally  equivalent  to 
0.23  •  gSL  for  each  time  period  of  the  walks  where  the  gSL  is  the  solar 
load  (W/m2)  as  generally  reported  in  the  weather  field  stations. 
Therefore,  the  net  solar  load  on  the  individual  is  optimally  expressed 
by  either  Rsoi  or  ERF  (W/m2),  and  both  functions  can  be  quantified 
during  transients  in  algebraically  different,  but  equivalent  heat  transfer 
methods,  as  described  below. 

During  field  operations,  Rsoi  (W/m2)  can  be  calculated  from 
0.835  ■  ERF.  If  Rsoi  (in  W/m2)  is  known  by  direct  pyranometer  mea¬ 
surements,  MRT  (°C)  can  be  estimated  from  0.0965 -(Rsoi  Ad)  +  Ta 
or  vice  versa  if  the  MRT  is  recorded  from  a  black  or  khaki-colored 
6-in.  (15.2  cm)  globe  or  4tt  radiometer  (15,  16,  25).  When  rough- 
estimated  weather  station  gSL  (W/m2)  values  are  the  only  available 
sources,  ERF  can  be  predicted  from  0.23 -gSL,  and  MRT  can  be 
determined  from  the  calculated  ERF  as: 

MRT(°C)  =  Ta  +  ERF/(feff-hr)  (9) 

Equation  9  therefore  sums  the  effect  of  mixed  solar  radiation  from 
direct,  diffused,  and  ground  radiant  loads  during  a  given  transient 
solar  radiation.  The  equation  can  also  be  used  to  integrate  radiant 
exchange  with  the  body  surface  when  MRT  differs  with  Ta  and 
estimate  the  effect  of  heat  storage  (S)  as  the  individual  exercises  over 
a  given  period  of  time  throughout  variable  solar  radiation  intensi¬ 
ties  (15). 

Heat  transfer  analyses.  In  the  present  study,  clothing  heat  and 
evaporative  potential  parameters  were  determined  using  a  regionally 
heated,  articulated  manikin  at  various  wind  speeds.  Each  element  of 
the  comprehensive  heat  balance  equation  was  analyzed  from  the  raw 
data,  specific  clothing  factors  by  applying  separate  y-coefficients  for 
thermal  and  vapor  resistances  in  an  articulated,  walking  manikin  (16, 
17,  24). 

The  raw  data  were  put  together  in  a  unified  spreadsheet  (Microsoft 
Excel)  for  later  analysis.  The  techniques  to  estimate  the  heat  transfer 
equations  above  and  all  other  analytical  equations  formulated  in  a 
previous  report  were  also  applied  to  this  study  (7,  17). 

Statistical  analyses.  Data  were  analyzed  using  a  variety  of  statis¬ 
tical  modules  (Statistica,  version  10,  Tulsa,  OK  and  STATA,  College 
Station,  TX),  the  latter  for  testing  and  programming  specific  comput¬ 
erized  bootstrap  analyses  and  Monte  Carlo  approaches  (12,  33,  34,  40, 
44,  45).  Simple  or  multiple  stepwise  regressions  analyzed  the  depen¬ 
dence  between  variables.  Since  some  subjects  in  night  trials  1-5  were 
used  in  repeated  studies,  the  differences  between  means  were  also 
analyzed  using  paired  f-tests  to  check  intraindividual  sweat  loss 
observations.  For  more  than  two  tests,  we  used  a  one-way  ANOVA 


for  repeated  measures  (26).  Otherwise,  the  values  were  tested  by 
unpaired  f-tests  over  the  various  trials.  The  assumptions  of  normality 
and  homogeneity  of  variance  for  parametric  procedures  were  checked 
using  Kolmogorov-Smimov  test.  When  the  assumptions  of  normality 
or  homogeneity  of  variance  were  not  met,  we  used  equivalent  non- 
parametric  tests  (Rruskal- Wallis  tests)  for  comparison  of  means  of 
two  independent  samples  (26).  Following  ANOVA,  the  data  were 
tested  using  Tukey’s  post  hoc  for  honestly  significant  differences  for 
unpaired  data  or  by  applying  a  Bonferroni  post  hoc  analysis  (26). 

A  composite  analysis  was  performed  using  pooled  sweat  rate  data 
from  77  subjects  completing  the  night  marches  and  550  laboratory 
observations  from  a  previous  study  (17).  The  data  were  parsed  and 
aggregated  into  a  separate  bootstrapping  analysis  to  determine  resid¬ 
ual  variations  distinct  in  the  original  data  (45).  In  the  bootstrap 
analysis,  the  hypothesis  is  that  a  given  sample  S  =  {Xj,  X2,  ■  ■  ■  ,X„ } 
(i.e.,  77  subjects  from  night  runs)  is  a  unique  cohort  of  the  given 
population  P  =  {xi,  X2, . .  .  ,Xn],  where  the  population  includes  all 
msw  data  from  the  original  PW  equation  analysis  (12,  17).  Efron  (12) 
developed  the  critical  statistic,  T*bs  =  t(S)  that  is  useful  in  analyses 
of  small  sample  bootstrapping.  This  statistic  supesedes  the  classic 
Student’s  t  to  test  small  departures  from  normality.  T*bs  was  deter¬ 
mined  as  an  estimate  of  the  corresponding  population  parameter  0  = 
t(P).  In  this  analysis,  0  is  assumed  as  a  vector  (bold  highlight)  of  all 
parameters  and  T*bs  the  corresponding  vector  of  estimates  (12,  45). 

The  bootstrap  analogy  using  Efron’s  criterion  ( 12)  in  our  prediction 
model  analysis  is  to  derive  T*bs  as  high  as  possible;  if  too  low  (<10 
or  so),  the  S  (individual  msw  output  value)  is  assumed  to  be  not  very 
well  predicted  by  the  purported  equation  or  contributes  too  large  an 
outlier  of  the  P  (or  too  far  out  of  the  domain  of  confidence  limits).  The 
bootstrap  estimate  of  the  standard  error  (SE)  of  T*bs  is  subsequently 
determined  as  the  square  root  of  the  bootstrap  variance  of  T*bs  as 
specified  by  Efron’s  criteria  (12). 

A  bootstrapped  correlation  matrix  was  next  created  using  output 
predictions  generated  from  Monte  Carlo  analysis  applying  the  prin¬ 
ciples  suggested  by  Picard  (33,  34).  The  Monte  Carlo  analysis  was 
followed  by  a  measure  of  goodness  of  fit  of  the  bootstrap  to  explain 
the  uncertainty  for  the  range  of  values  of  the  pooled  data.  The 
goodness  of  fit  test  that  was  used  in  our  study,  post  Monte  Carlo 
analysis  was  the  Schwarz-Bayesian  criteria  (SBC)  (40),  which  is  a 
function  of  the  natural  log  [ln(nx)]  of  the  number  of  observations  n, 
the  sum  of  squared  errors  (SSE),  and  the  number  of  independent 
variables  k^p  +  1  where  k  includes  the  intercept  as  shown  in  Eq.  10: 

(SBC)  =  n- In  (SSE/w)  +  k  In  n  (10) 

Equation  10  thus  comprises  the  summed  square  of  residuals  SSE  = 
Sum(,=1  to  „)[w,-(y,-  —  fd2],  where  y is  the  observed  msw  data  value  and 
f  is  the  predicted  (PW)  value  from  the  fit;  and  w,-  is  the  weighting 
applied  to  each  data  point,  typically  w,  =  1. 

The  bootstrapping  analysis  enabled  determination  of  an  estimated 
standard  error  (SE)  of  correlation  of  the  modified  PW  equation  (12), 
as  noted  above,  using  measured  sweat  rates  in  the  new  night  field 
cohort  data,  assuming  such  data  sets  came  from  an  anomalous  mul¬ 
tivariate  population  (12,  34)  (i.e.,  uncertainty  due  to  dissimilar  field 
data  tested  vs.  conventional  indoor  chamber  data).  Individual  sweat 
rate  values  from  the  operational  marches  (trials  1-7)  were  analyzed 
using  the  data  splitting  technique  (33).  The  total  data  set  comprised  all 
observed  sweat  rate  values  from  individuals  completing  each  trial,  but 
trials  1-5  (Table  1,  night  marches)  were  analyzed  independently  from 
trials  6  and  7  (Table  2,  transient  solar  marches).  All  data  are  expressed 
as  means  ±  SD  or  as  means  ±  95%  prediction  intervals. 

RESULTS 

Night  trials  1-5.  Table  4  provides  subject  anthropometric 
characteristics,  march  distance,  load,  and  selected  physiologic 
responses  for  the  seven  trials.  During  trials  4  and  5,  the 
calculated  oxygen  uptake,  mean  weighted  skin  temperature, 
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Table  4.  Key  descriptive,  environmental,  and  physiological  data  from  field  experiments  used  in  current  validation 


Trial 

Ta,  °C 

Pw,  kPa 

V,  m/s 

Distance,  km 

Load,  kg;  L/mass,  % 

BSA,  m2 

Body  Mass,  kg 

WC 

V02,  1/min 

SR,  1/h 

i 

25.1 

2.28 

0.6 

24.0 

30.9  ±  6.7;  40.4  ±  8.6% 

1.94  ±  0.09 

76.7  ±  5.5 

33.81  ±  0.26 

2.42  ±  0.21 

1.02  ±  0.21 

2 

20.8 

1.95 

0.5 

29.0 

30.4  ±  7.2;  40.2  ±  10.7% 

1.94  ±  0.13 

77.5  ±  9.7 

31.29  ±  0.41 

2.43  ±0.18 

0.97  ±  0.22 

3 

17.4 

1.30 

1.3 

37.0 

35.2  ±  5.9;  45.8  ±  8.1% 

1.94  ±  0.13 

77.4  ±  9.5 

29.17  ±  0.67 

2.72  ±  0.30 

0.72  ±  0.23 

4 

17.4 

1.17 

0.9 

14.0 

15.1  ±  4.5;  21.5  ±  5.1% 

1.85  ±  0.10 

69.9  ±  6.2 

29.76  ±  0.93 

1.84  ±  0.21 

0.57  ±  0.16 

5 

20.0 

1.65 

1.1 

14.0 

17.0  ±  3.3;  24.4  ±  3.5% 

1.85  ±  0.10 

69.7  ±  6.3 

30.91  ±  0.58 

2.16  ±  0.21 

0.68  ±  0.18 

6 

25.9 

1.81 

1.3 

16.4 

15.3  ±4.3;  21.1  ±6.3% 

1.89  ±  0.18 

73.7  ±  10.8 

33.11  ±  0.33 

1.69  ±  0.23 

0.99  ±  0.17 

7 

16.2 

0.94 

1.4 

18.3 

21.0  ±  8.2;  30.8  ±  11.3% 

1.83  ±  0.10 

67.9  ±  6.4 

29.30  ±  0.96 

1.68  ±  0.22 

0.61  ±  0.14 

All  data  are  means  ±  SD.  Air  temperature  (Ta),  ambient  air  movement  (V),  ambient  water  vapor  pressure  (Pw),  distance  traversed  (km),  load  carried  (L,  kg), 
body  surface  area  (BSA),  body  mass  (kg),  mean  skin  temperature  (TSk),  and  observed  sweating  rate  (SR)  are  end  values. 


and  absolute  sweat  loss  values  were  lower  ( P  <  0.01)  com¬ 
pared  with  mean  values  in  trials  1  and  2  (Table  4).  The  loads 
carried  during  trials  1,  2,  and  3  were  larger  ( P  <  0.01)  than  the 
other  trials.  The  load  (L)/body  mass  (L/kg  X  100)  percentage 
carried  in  trials  1,  2,  and  3  (ranges  40-46%,  P  <  0.05)  was 
higher  than  %L/body  mass  apparent  in  the  other  two  night 
trials. 

Table  5  provides  measured  sweat  rates  for  the  night  trials 
1-5,  the  respective  PW  output  (17),  and  the  Shapiro  (41) 
predicted  msw  values.  The  average  deviation  of  each  prediction 
output  compared  with  measured  sweat  rate  is  provided  (A  of 
model  —  measured).  For  all  night  trials,  the  deviation  from  PW 
output  values  was  not  different  from  the  measured  sweat  rates 
(mean  ±  SD  =  28  ±  9.6  g-m_2-h  1  for  trials  1-5).  For  all 
night  trials,  the  Shapiro  (41)  values  deviation  was  greater  ( P  < 
0.01  to  P  <  0.001)  than  the  measured  sweat  rates. 

Figure  1  presents  the  bootstrap  analysis  of  the  PW  output 
values  (as  dependent  variable)  plotted  against  measured  sweat 
rates  for  trials  1-5  (77  experiments)  and  the  previously  col¬ 
lected  laboratory  data  (17).  The  analysis  employed  analytic 
procedures  using  the  bootstrap  analysis  on  matched  residual 
variances  as  described  in  methods  (33,  45).  This  analysis 
incorporates  Ereq  and  Emax  parameters  implemented  from  PW 
prediction  of  msw.  Thus  the  present  study  data  served  for  both 
prediction  and  cross-validation  of  the  original  fitted  PW  equa¬ 
tion  (17,  33). 

The  bootstrap  parameter  regression  of  this  analysis  is  shown 
in  Fig.  1.  The  sweat  rates  calculated  from  the  heat  exchange 
variables  (Ereq  and  Emax)  of  trials  1-5  aggregated  (closed 
triangles)  with  the  PW  equation  vs.  the  sweat  rates  from  the 
ARIEM  lab  data  (open  circles)  had  outliers  that  predicted 
lower  than  the  95%  prediction  line  in  9  cases  (=1.5%  of  over 
550  data  points).  Conversely,  there  were  12  cases  of  outliers 
predicting  msw  higher  than  95%  confidence  limit  (=2%  of  the 
total  data  points).  The  higher  predicted  sweat  rates  are  within 
the  domain  of  variability  of  the  bootstrapping  variances  inher¬ 
ent  in  the  PW  algorithm  (12,  45)  and  the  whole  regression 


equation  displayed  a  SEE  =  ±69.2  g-m_2-h_1.  After  1,000 
iterations,  the  bootstrap  correlation  (rbS)  parameter  estimate 
defaulted  to  0.81  with  a  standard  error  of  the  |3-correlation 
coefficient  of  ±0.014  [Efron’s  (12)  T*bs-value  =  31.3,  P  < 
0.0001]  for  the  pooled  data.  The  data  now  include  field  night 
run  data  that  were  combined  with  the  original  laboratory 
experiments.  The  sweat  rate  data  exhibited  essentially  no 
heteroscedasticity  (unequal  variances  among  the  groups)  (26, 
33,  34,  45).  Following  the  bootstrapping  analysis,  a  Monte 
Carlo  parsing  (34)  was  next  executed  on  the  combined  data  set 
to  determine  the  optimal  model  selection  based  on  an  uncer¬ 
tainty  goodness  of  fit  of  the  PW  equation  using  SBC.  For  the 
final  Monte  Carlo  iterations,  a  SBC  of  0.725  was  derived  (1  = 
perfect  fit  to  the  data)  that  suggests  that  the  PW  model  using  all 
combined  data  values  (night  +  laboratory  studies)  still  exhibits 
small  residual  random  error  components  (0.275).  As  previously 
mentioned,  a  greater  proportion  of  variance  causing  over-  or 
underprediction  errors  is  still  accounted  for  in  the  lumped 
proportional  control  (^-coefficients)  parameters  present  in  the 
original  fuzzy  piecewise  (PW)  model  (17).  The  data  fit  also 
conformed  with  the  less  robust  Ordinary  Least  Squares  (OLS) 
(r  =  0.71,  P  <  0.001)  and  confirms  that  sweat  loss  prediction 
using  the  PW  model  was  not  adversely  limited  by  any  outliers 
changing  the  gain  algorithm  ((3-coefficients)  in  the  model 
prediction  structure  due  to  too  high  or  too  low  sweat  rates 
among  the  subjects  during  the  various  trials  1-5  (26,  33,  46). 
We  next  analyzed  sweat  rates  in  the  marches  where  solar  load 
was  apparent. 

Transient  solar  radiation  trials  6-7.  Table  4  displays  the 
subject  anthropometric  characteristics,  load,  and  selected  phys¬ 
iologic  responses  for  trials  6  and  7.  In  trial  6,  20  subjects  began 
marching  at  0330,  with  explicit  solar  gain  appearing  at  0500 
and  terminating  with  a  maximum  global  solar  load  (gSL)  of 
341  W/m2.  During  this  march,  Ta  was  25.9°C  but  the  effective 
heat  stress  gain  owing  to  the  increased  MRT  was  35.9°C.  As 
shown  in  Table  2,  trial  6  was  a  scenario  that  comprised  two 
marching  tasks  up  a  mountainous  gravel  road  (median  terrain 


Table  5.  Measured  sweat  rate  compared  with  predicted  output  during  night  marches  of  trials  1-5 


Trial  1 

Trial  2 

Trial  3 

Trial  4 

Trial  5 

No.  of  subjects 

13 

12 

12 

20 

20 

SR,  g-m_2-h_1 

524  ±  96 

502  ±  121 

366  ±  104 

307  ±  84 

367  ±  96 

PW  prediction  (Ref.  17) 

493  ±  75 

489  ±  55 

405  ±  41 

338  ±  28 

393  ±  60 

Shapiro  prediction  (Ref.  41) 

854  ±  113 

847  ±  106 

752  ±  42 

636  ±  30 

727  ±  101 

A(PW  -  SR) 

31 

13 

39 

31 

26 

A(Shapiro  —  SR) 

330* 

345* 

386f 

329* 

360t 

Values  are  means  ±SD.  Comparison  of  average  deviation  A  of  model  output  to  actual  sweat  rate  (SR):  *P  <  0.01;  tP  <  0.001;  A(PW  —  SR)  in  all  other 
trials,  not  significant. 
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Sweat  Rate,  msw  (g-rrf2-h'1) 

Fig.  1 .  PW  prediction  plotted  as  a  function  of  sweat  rate  during  night  trials  1-5 
(closed  triangles)  aggregated  with  previous  laboratory  results  (open  circles) 
(17).  Data  include  bootstrapping  correlation,  tbs  (12,  34,  45);  solid  line  is  the 
identity  line,  PW  =  sweat  rate.  The  bootstrap  estimate  of  standard  error  (SE) 
for  the  pooled  data  was  ±0.014  with  a  T*b  =  31  (12).  The  standard  error 
estimate  of  the  regression  (SEE)  was  ±  69.2  g-m_2-h_1.  Typically,  once  the 
original  fuzzy  piecewise  PW  equation  is  applied  (17,  44),  this  output  value  can 
be  substituted  in  the  above  OLS  equation  (125  +  0.67 -PW)  to  calculate  the 
msw  required;  for  example,  600  g-m_2-h-1  determined  from  the  equation 
yields  522  g-m_2-h_1  msw,  within  the  95%  prediction  interval  of  the  boot¬ 
strap.  These  data  now  include  field  night  runs  combined  with  the  original  lab 
experiments. 

coefficient,  t|  =  1.3)  with  work/rest  cycles  of  65  min  work  and 
18  min  rest  followed  by  another  march  of  48  min  work/ 15  min, 
and  concluding  with  a  resting  period  of  62  min.  Total  exposure 
time  to  the  transient  solar  load  was  208  min. 

In  trial  7,  13  subjects  marched  on  a  mixed  paved  and  gravel 
road  (t|  =  1.25).  These  marches  began  at  1420  at  which  time 
maximum  global  solar  load  (gSL)  was  612  W/m2.  During  this 
march,  Ta  was  16.2°C  and  the  resulting  MRT  was  26.2°C.  The 
work/rest  marches  in  trial  7  continued  until  1 800  at  which  time 
solar  load  reached  a  nadir  with  minimal  or  no  solar  load.  Trial 
7  comprised  220  min  of  exposure  to  a  transient  solar  load 
which  had  two  separate  marches  with  work/rest  cycles  that 
included  58  min  work/15  min  rest  and  another  session  that  was 
71  min  work/14  min  rest,  followed  by  a  final  resting  phase  of 
62  min. 

Figure  2  shows  the  mean  (±SD)  solar  radiation  (W/m2)  as  a 
function  of  time  of  day  during  the  marches  using  Rsoi  derived 
from  gSL  for  trials  6  and  7.  The  MRT  calculated  from  Rsoi  and 
ERF  increased  from  a  value  of  25.7°C  to  as  high  a  50°C  during 
trial  6,  and  decreased  from  58°C  to  as  low  as  16.1°C  during 
trial  7.  Rsoi  and  ERF  as  explained  in  methods  incorporated 
each  timed  output  to  determine  Eieq  and  Emax  for  later  input  in 
the  modified  PW  prediction  equation  in  trials  6  and  7. 

Based  on  Ereq  and  Emax  adjustments  using  Rsoi,  MRT,  and 
Veff  calculated  for  the  transient  solar-night  trials,  a  revised  PW 
(Pw,sol)  ordinary  least  squares  (OLS)  equation  for  transient 
solar  radiation  conditions  was  developed  [msw  =  0.72- 
(/Tv, sol)  +  135],  where  Pw,sol  reflects  the  original  PW  equa¬ 
tion  (17)  with  Ereq  now  modified  to  include  Rsoi  (Eieq.Soi).  The 
value  of  Ereq,soi  must  be  first  determined  by  pertinent  solar 
coefficients.  Table  6  provides  the  results  of  a  multiple  regres¬ 
sion  analysis  conducted  using  sweat  rate  (dependent  variable) 
from  each  of  subjects  completing  trial  6  and  7.  The  output  of 


three  prediction  equations  were  compared:  1 )  Pw,sol  as  devel¬ 
oped  in  the  current  analysis;  2)  Shapiro  (42),  which  has  been 
discussed;  and  3)  HSDA,  which  is  often  used  by  the  military  to 
predict  sweat  rate  for  outdoor  activities  (3,  31).  The  magnitude 
of  the  |3-coefficients  for  each  prediction  equation  describes  the 
relative  contribution  of  each  independent  variable  in  the  pre¬ 
diction  of  sweat  rate.  Table  6  indicates  that  Pw,sol  had  the 
highest  (3-coefficients  and  was  the  only  equation  that  accounted 
for  significant  (P  <  0.0001)  variance  contributions  from  each 
predictor  variable  in  the  orthogonal  correlation  structure  of  the 
multiple  regression  analysis  (26).  The  Pw,sol  output  was  also 
significantly  correlated  with  measured  sweat  rates  from  the 
pooled  transient  solar  experiments.  Figure  3  also  presents 
output  showing  the  OLS  analysis  of  the  predictability  of 
/Tv,  so  I  output  values  (here  plotted  as  the  independent  variable) 
and  measured  msw  for  trials  6  and  7.  A  significant  correlation 
(r  =  0.96;  P  <  0.001)  and  SEE  (±33.8  g-m_2-h_I)  were 
detected.  Once  the  equation  is  solved  for  the  various  indepen¬ 
dent  variables  in  the  heat  exchange,  the  model  output  appears 
as  a  simple,  practical  tool  in  prediction  of  msw  for  a  given 
transient  solar  load. 

Additional  data  exploration.  Further  investigation  of  the 
combined  data  allowed  modification  of  the  laboratory-derived 


Fig.  2.  Rsoi  output  values  (25)  plotted  vs.  time  of  day  during  the  transient  solar 
experiments  of  trial  6  (top  panel)  and  trial  7  (bottom  panel).  Saw-toothed  dips 
in  Rsoi  of  trial  6  are  due  to  mixed  cloud  formation  impacting  the  solar  flux  in 
some  trials. 
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Table  6.  Multiple  regression  of  sweat  rate  (g-m~2 -h~2)  vs.  output  from  Shapiro  equation  ( Ref.  42)  (Shap_out),  Heat  Strain 
Decision  Aid  (Refs.  3,  31)  (HSDA),  and  PW,sol  equation  for  pooled  data  from  the  transient  solar  trials  6  and  7 


Model  Equation 

p* 

SE  of  p* 

P 

SE  of  |3 

r(23) 

P  < 

HSDA 

-0.020 

0.062 

-0.113 

0.350 

-0.32 

0.748 

Pw,sol 

0.929 

0.066 

0.66 

0.047 

14.05 

0.0001 

Shap_out  (Ref.  42) 

0.066 

0.075 

0.21 

0.240 

0.87 

0.391 

|3  *  is  standardized  regression  coefficient;  [i  is  raw  regression  coefficient;  magnitude  of  (j  coefficients  compares  relative  contribution  of  each  independent 
variable  in  prediction  of  the  sweat  rate  for  each  subject.  HSDA  out  (Refs.  3,  31)  and  Shap_out  (Ref.  42)  equations  were  not  good  predictors  of  sweat  rate  (P  > 
0.05,  NS)  for  the  transient  solar  experiments.  N  =  33  subjects. 


PW  equation  with  a  focus  on  cool  environments  during  resting 
or  low  metabolic  activities.  Cool  (<16-18°C)  ambient  tem¬ 
peratures  result  in  low  skin  wettedness,  («  <  0.15)  with  water 
loss  primarily  from  skin  diffusion  (5,  15,  19,  22).  This  effect 
takes  place  often  when  Emax  is  >400  g  •  m  2  ■  h  1  (approxi¬ 
mately  >270  W/m2)  and/or  Ereq  becomes  <50  g-m~2-h_1 
(13,  20,  22,  41).  In  such  environments,  the  equation  to  predict 
sweat  rates  (within  131  g  •  m  2  •  h  1  of  actual  values)  can  now 
be  determined  from  the  following: 

PW,  cool(g  •  m”2  •  h  ')  =  350  +  (U) 

1.527  •  Ereq  -  0.87  •  Emax 

The  above  equation  does  not  alter  the  (3-coefficients  (algorithm 
gain)  from  the  original  PW  equation  but  increases  the  intercept 
by  203  g  ■  m  ~2  ■  h  1  when  Ereq  is  low,  but  Emax  is  high.  There¬ 
fore,  for  the  Ereq  and  Emax  mentioned,  Eq.  10  predicts  msw  = 
~2.5  g/min  occurring  with  elevated  wind  speeds  during  labo¬ 
ratory  and  night  (limited  outdoor  field)  experiments. 

DISCUSSION 

This  study  validated  a  previously  developed  sweat  rate 
prediction  equation  and  derived  pertinent  biophysical-physio- 
logical  corrections  to  modify  that  equation  for  the  valid  pre¬ 
diction  of  sweat  rate  (water  needs)  during  transient  solar 
radiation  and  cool  conditions.  The  new  Pw,sol  and  PW,cool 
equations  have  important  implications  for  military  and  disaster 
relief  planning  and  a  variety  of  other  occupational  problems 
related  to  the  logistical  planning  for  water  needs.  We  also 


Trials  6  &  7 


Pw,sol  (g-m'V) 

Fig.  3.  Sweat  rate  from  trials  6  and  7  plotted  as  a  function  of  Pw,sol  equation 
derived  from  heat  exchange  parameters  following  analytical  solution  of  Ereq,  soi 
as  explained  in  text. 


demonstrated  that  several  other  often  used  sweat  prediction 
equations,  for  outside  activities  with  a  solar  radiation,  were  far 
less  accurate  compared  with  the  modified  /V.sol  equation. 

During  the  transient  solar  experiments,  we  observed  high 
sweat  rates,  despite  the  low  ambient  temperatures,  due  to  the 
impact  of  the  solar  gain  quantified  by  elevated  MRT.  Even 
though  load  carriage  and  metabolic  level  would  be  considered 
nominal,  the  ambient  temperatures  (Ta  =  25.9°C  in  trial  6  and 
16.2°C  in  trial  7)  were  influenced  by  the  added  rise  in  MRT 
that  exacerbated  the  sweat  rates  found  in  trial  6  compared  with 
the  values  in  trial  7.  During  the  transient  solar  gain,  average 
sweat  rates  calculated  from  the  20  individuals  in  trial  6  were 
523.4  gm~2-h-1  (±73.6  SD),  compared  with  332.5 
g-nW2-h-1  (±69.4  SD)  in  the  13  individuals  of  trial  7  (58% 
increase,  P  <  0.001,  nonpaired  t  =  7.45).  These  values  amount 
to  some  16.6  g/min  higher  sweat  rates  observed  in  trial  6  vs. 
10.5  g/min  in  trial  7.  The  necessary  water  requirements  to 
sustain  such  sweat  rates  to  maintain  thermoregulation  are 
greatly  increased  by  the  solar  impact  (ERF  and  MRT)  during 
the  walks  despite  the  mild  ambient  conditions.  Interestingly  a 
constant  cool  ambient  temperature  (as  in  trial  7;  Table  4), 
reducing  solar  radiation  from  a  maximum  (612  W/m2)  in  the 
late  afternoon  to  zero  at  night,  markedly  reduced  sweat  rates 
(and  water  needs,  respectively)  compared  with  early  morning 
marches  in  which  solar  load  gradually  accumulates  with  a 
moderate  air  temperature.  These  findings  are  consistent  with 
Gagge  and  Hardy  (14)  who  first  observed  that  thermoregula¬ 
tory  sweating  (evaporative  heat  loss)  changes  were  highly 
correlated  with  ERF  altered  by  short  bursts  of  infrared  heating 
(heat  flux  from  IR  quartz  lamps)  at  Ta  =  23-24°C.  They  also 
proved  that  as  ERF  increased  monotonically  at  each  higher 
constant  ambient  temperature  (from  9°C  to  30°C),  the  evapo¬ 
rative  heat  loss  by  sweating  increased  proportionately.  Nielsen 
et  al.  (30)  also  showed  similar  high  sweat  rates  during  outdoor 
cycle  exercise  in  the  sun. 

The  present  study  confirmed  our  initial  hypothesis  that  sweat 
rate  predictive  equations  developed  from  laboratory  studies 
(17)  are  suitable  for  matching  sweat  requirements  (or  at  least 
were  not  statistically  overpredictive;  i.e.,  SEE  <  ±70 
g  ■  m  " 2  ■  h  ~ 1 )  to  actual  sweat  loss  for  disparate  military  activi¬ 
ties  during  the  night  time.  This  was  an  anticipated  finding  since 
the  experiments  involved  similar  conditions  (Ta,  RH,  clothing 
thermal  resistance  values,  and  work  rates)  as  in  previous 
environmental  chamber  studies  (7,  17).  Additionally,  the  major 
finding  from  the  present  study  was  improvement  of  a  sweat 
loss  prediction  equation  for  transient  solar  loads  quantified  by 
Rsoi,  ERF,  and  MRT.  To  our  knowledge,  this  has  not  been 
possible  before  and  may  require  additional  investigation  par¬ 
ticularly  at  higher  ambient  (hot  dry  and  tropical)  environments. 
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Most  experiments  that  involve  solar  load  challenges  in  humans 
are  done  during  constant  thermal  radiation  exposures  and  brief 
time  periods  (14,  30).  In  that  respect,  the  avenues  of  heat 
exchange  (Eieq  and  Emax)  are  assumed  as  being  under  steady- 
state  conditions.  Our  data  analysis  shows  that  Rsoi  as  a  variable 
input  into  the  heat  balance  equations  can  be  used  for  estimating 
sweating  during  transient  solar  environments. 

The  prediction  equations  for  estimating  sweat  rates  using 
Rsoi  garnered  from  the  present  experiments  are  an  important 
new  finding  because  they  include  persistent  effects  of  heat 
storage  (S)  as  core  and  skin  temperatures  increase  (37)  and 
thereby  integrate  complex  solar  environments  more  rigorously 
than  previous  investigations  (4,  42).  We  also  derived  by  ex¬ 
ploration  of  present  and  past  (17,  41)  data  sets,  an  adjusted  PW 
equation  (PW,cool)  that  is  applicable  for  cooler  weather  con¬ 
ditions  when  PW  might  predict  negative  sweating  values  be¬ 
cause  Ereq  is  too  low  or  Emax  is  too  high.  The  adjustment 
requires  only  a  change  to  the  legacy  PW  equation  intercept 
term.  We  recommend  that  these  new  equations  should  be 
included  in  current  thermal  models  (3,  24,  46)  especially  for 
field  use  when  solar  load  is  present.  Fundamentally,  once  the 
various  discrete  parameters  (Veff,  Rsoi,  Wclo,  etc.)  influencing 
dynamic  clothing  heat  transfer  and  walking  factors  are  quan¬ 
tified  into  the  heat  balance  determining  Ereq  and  Emax,  the 
Pw,sol  equation  is  valid  in  predicting  msw  (within  an  SEE  ±  34 
g  ■  m  ‘ 2  •  IT  1  or  1.1  g/min  of  actual  sweat  rate  values  or  predict¬ 
ing  fluid  needs  =  ±0.1  liter/h)  (Fig.  3).  The  approximation  is 
a  best  case  prediction  scenario  possible  in  field  studies  that  is 
limited  by  the  domain  of  thermal  environments,  metabolic 
intensities,  and  clothing  ensembles  studied  here  considering 
the  sweat  rate  variability  of  a  general  population  (18). 

Water  needs  are  currently  estimated  for  the  US  military 
using  HSDA  or  similar  approaches  (3,  21,  23,  28).  Our  past  (7, 
17)  and  this  study  clearly  demonstrate  that  the  PW  equations 
provide  for  more  accurate  predictions  of  sweat  rates  during 
various  exercise  intensities  and  environments.  A  caveat  of  the 
HSDA  guideline  is  that  the  input  value  for  skin  temperature 
executed  in  the  calculation  of  dry  heat  loss  (R  +  C)  is 
determined  by  a  fixed  value  of  36.5°C  [skin  saturation  vapor 
pressure  of  45  Torr  (5.99  kPa)]  rather  than  actual  values  as 
used  in  this  study.  Furthermore,  the  HSSk  values  employed  that 
determine  MRT  in  that  guideline  are  highly  sensitive  to  fluc¬ 
tuations  (y-values)  occurring  in  a  variety  of  thermal  and  vapor 
resistances  (20,  22,  24)  as  determined  using  thermal  manikins. 
For  example,  Table  6  showed  that  calculation  of  the  HSSk  factor 
in  the  HSDA  model  (3)  did  not  predict  sweat  rates  reliably 
during  the  transient  solar  trials  in  this  study.  One  reason  may 
be  because  Hssk  is  a  collective  equation  in  the  HSDA  empirical 
model  that  incorporates  MRT,  but  is  coupled  to  arbitrary 
clothing  and  air  movement  coefficients  (3).  Therefore,  any 
derived  sweat  loss  prediction  can  therefore  be  too  high  or  too 
low  based  on  the  inaccuracies  of  the  Hssk  factor.  A  new 
rational  multimodel  also  demonstrated  large  residuals  even 
when  applied  to  group  data  (46).  The  variety  of  equations:  PW, 
Pw,sol,  and  PW,cool  could  be  implemented  in  various  model 
input  interfaces  (for  example  by  using  a  Boolean  logic  “do- 
loop”  contingency)  when  Ta  is  less  than  20°C,  solar  load  is 
present,  wind  movement  is  high,  or  cloudy  outdoor  activity  is 
planned  (provided  all  Ereq/Emax  caveats  and  domain  of  validity 
are  upheld). 


In  summary,  we  performed  a  series  of  outdoor  military  field 
studies  in  which  sweat  rate  was  quantified  in  110  soldiers 
exercising  during  night  trials  and  during  transient  solar  radia¬ 
tion  conditions.  We  validated  our  recently  published  sweat  rate 
prediction  equation  (17)  for  outdoor  military  activities  and 
modified  it  so  it  predicts  msw  accurately  for  transient  solar 
radiation  and  cooler  ambient  conditions.  We  compared  our  PW 
sweat  prediction  equations  to  other  existing  predictions  equa¬ 
tions  and  established  that  PW  was  an  improved  calculation 
readily  applicable  during  field  training  exercises  with  (Pw,sol), 
without  (original  PW)  transient  solar  loads,  and  finally  when 
environmental  conditions  embrace  low  Ereq  and  high  Emax 
capacities  (PW,cool)  during  moderate  exercise  intensities.  The 
family  of  PW  sweat  prediction  equations  allows  easy  imple¬ 
mentation  tools  in  civilian,  sports,  and  military  medicine  com¬ 
munities  needed  to  estimate  water  needs  during  a  variety  of 
heat  stress/exercise  conditions. 
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